Cosmologies in Horndeski's second-order 
vector-tensor theory 



o 



John D. Barrow," Mikjel Thorsrud fe and Kei Yamamoto ab 

a DAMTP, University of Cambridge, 

Wilberforce Road, Cambridge, CBS OWA, United Kingdom 
b Institute of Theoretical Astrophysics, University of Oslo, 

P.O. Box 1029 Blindern, N-0315 Oslo, Norway 

E-mail: J . D . Barrow@damtp . cam . ac . uk, mikj el . thorsrud@astro . uio . no, 
K . Yamamoto@damtp . cam .ac.uk 

Abstract: Horndeski derived a most general vector-tensor theory in which the vector 
field respects the gauge symmetry and the resulting dynamical equations are of second 
order. The action contains only one free parameter, A, that determines the strength of the 
non-minimal coupling between the gauge field and gravity. We investigate the cosmological 
consequences of this action and discuss observational constraints. For A < we identify 
singularities where the deceleration parameter diverges within a finite proper time. This 
effectively rules out any sensible cosmological application of the theory for a negative non- 
minimal coupling. We also find a range of parameter that gives a viable cosmology and 
study the phenomenology for this case. Observational constraints on the value of the 
coupling are rather weak since the interaction is higher-order in space-time curvature. 
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1 Introduction 

In the past few decades, modifying the Einstein's theory of gravitation has been an active 
area of research [1], driven chiefly by the search for different varieties of inflation, the desire 
of some to explain flat galaxy rotation curves without dark matter [2], and the challenge of 
explaining why the expansion of the universe started to accelerate at late times. It is also 
natural to question the validity of general relativity, not least because of its ultra-violet 
behaviour which does not give a well-defined quantum field theory. There is a growing 
prospect of testing any such deviations from general relativity in very strong gravity fields 
by searching for the signatures of gravitational waves created by high-energy astrophysical 
phenomena, such as black hole mergers, or by scrutinising detailed observations of the mi- 
crowave background anisotropy and statistics in the light of particular theories of inflation. 
While there are a plethora of inflationary models, it is difficult to modify general relativity 
without spoiling its appealing features and typical modifications end up introducing new 
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scalar degrees of freedom, as was the case in the case of Brans-Dicke gravity [3] and its 
scalar-tensor generalisations [4-6], or f(R) lagrangian theories of gravity [7—11]. 

Maxwell's classical theory of electromagnetism has also been extremely well tested and 
there are already strict constraints on potential modifications such as a non-zero photon 
mass [12-15] or varying fine structure constant [16-24]. However, there has been renewed 
interest in modified electromagnetism in cosmology and there have been attempts to in- 
corporate its effects into the dynamics of early universe, particularly during inflation, or 
to provide an explanation for cosmological magnetic fields [25-40]. Simple vector fields 
themselves are known to have difficulties in producing inflation. The time variation of the 
vector field is governed by the covariant time derivative of the field. Since the Christoffel 
symbols for an expanding cosmological model are of order the Hubble expansion rate it 
is not possible for the vector field to satisfy a slow-roll condition in the way that a scalar 
field can [41]. However Einstein-aether theories offer an alternative that permits inflation 
[1, 24, 42-51]. In addition, although it was noted that a non-minimal coupling to the 
space-time curvature could drive accelerated expansion [52-55], these scenarios suffer from 
various instabilities created by additional degrees of freedom arising from the lost gauge 
symmetry [56-58]. A different type of extension of the Maxwell case is provided by the 
extension to Yang-Mills fields where there can be chaotic behaviour and arbitrarily low 
levels of anisotropy [59-61]. 

Recently, it has been mentioned that a vector-tensor theory, first proposed by Horn- 
deski in 1976 [62], could lead to an instability of conventional inflationary universe through 
the non-minimal coupling between the vector field and gravity [63]. The action was de- 
rived by demanding second-order dynamical equations that reduce to Maxwell's equations 
when evaluated on a Minkowski background and conservation of the U(l) current. These 
requirements result in only one additional term in the Lagrangian, and therefore a single 
free coupling parameter. Later, it was noticed that this theory falls into a special class of 
Kaluza-Klein reductions from higher-dimensional Lovelock invariants [64-66]. In contrast 
to the Horndeski scalar field theory [67], which has been discussed in attempts to construct 
the most general viable scalar-tensor theory recently, [68-73], except for a brief examina- 
tion of the static electromagnetism arising from this action [74], it appears to have escaped 
attention. Apart from being briefly mentioned in [63], its cosmological consequences have 
not been studied. 

In this paper, we will investigate the simplest cosmological model, which is well under- 
stood in the minimally coupled case of a Maxwell electromagnetic field [75, 76], containing a 
perfect fluid and a vector field whose dynamics are described by the Horndeski Lagrangian. 
We find the following results; 

1. The instability found in [63] for negative values of the coupling constant persists in 
the nonlinear regime and the universe eventually hits a singularity; 

2. For a positive coupling constant, the electric field can still be amplified during the 
radiation-dominated era while giving a viable cosmology subject to some constraints 
on the allowed expansion rate changes at the epoch of primordial nucleosynthesis. 
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The first result effectively rules out any interesting cosmological application of this 
theory with a negative coupling. For a positive coupling constant, the modification is 
rather tame and an enormous value of the coupling in units of the Planck mass is allowed 
because of the higher-order nature of the modified term. However, the dynamics is of a 
phenomenological interest. 

The article is organised as follows. In the next section, the theory is introduced and the 
modified Einstein-Maxwell equations are presented. Section 3 is the main part of the article 
where the dynamics of purely electric component are studied in an axisymmetric Bianchi 
type I universe. In section 4, we repeat the previous analysis for magnetic component. 
Section 5 discusses observational constraints and in section 6 we summarise our principal 
results. 

2 Horndeski's second-order vector-tensor theory 

In 1976, Horndeski showed that the general Lagrangian that can be constructed from a 
metric g ab and a vector field A a in four-dimensional space-time that satisfies the following 
conditions [62]: 

1. the field equations contain at most second-order derivatives of g ab and A a (and do 
contain a second-order term); 

2. the dynamical equations for A a respect charge conservation i.e. V a (dC/dA a ) = 0; 

3. the dynamical equations for A a reduce to Maxwell's equations when evaluated on 
Minkowski space-time; 

takes the following form: 

C = ^^g-R--^g~F ab F ab + C H , (2.1) 

where M p i is the reduced Planck mass, R is the Ricci scalar and F ab = d a A b — d b A a is 
the Faraday tensor. The last term is Horndeski's modification which can be expressed in 
several different ways as 

—gb ab e f kl F ab F^R cd kl (2.2) 

gF ab F cd * R * ab cd (2.3) 

—g (RF ab F ab - AR ab F ac F b c + R abcd F ab F cd ) . (2.4) 

The dimensionless non-minimal coupling constant A is the only parameter of the theory. 
Our aim is to investigate the cosmological consequences with an arbitrary value of A and 
to determine the parameter range yielding viable phenomenology. The other terms in the 
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Lagrangian (2.1) are normalized so that it reduces to the Einstein-Maxwell theory when 
A = 0. Using the Levi-Civita tensor r] abc d, we defined the generalised Kronecker's delta by 



rabcd rn srfr « rf ^ ^abcd^ 

6 efkl ~ o [e o f d k d q - - — rj rj e f k i 

and the double dual of Riemann by 

*Jf* ab — -n abe f n j, , F! kl 
* n * cd — Vcdkl-n-ef ■ 

Ch was later identified with the Lagrangian obtained by Kaluza-Klein reduction from the 
five-dimensional Lovelock invariant 

K = R abcd R abcd - AR ab R ab + R 2 

where R a bcd is the Riemann tensor in five dimensions [64]. 

In this paper we use the sign conventions of [77] for the metric, Ricci and Riemann 
tensors which are different from those adopted in the previous studies of this model [62, 74]. 
The dynamical equations derived from this Lagrangian are given as follows: 

Variation with respect to 

1 

4 

where 



Mlfiij = FiaFj a - 7 F ab F ab 9ij + Tij , (2.5) 



\ ( V a * F lb V b * F a 3 + F ab F c a * R* ibjc ) , (2.6) 



Variation with respect to Ai 



V a F ia - -^VaF bc * R* iabc = 0. (2.7) 



We define the dual Faraday tensor as usual: 

*F ab = ^VabcdF cd . 

In ref.[63], it was observed that (2.7) evaluated on a Friedmann-Lemaitre- Robertson- Walker 
(FLRW) background could lead to an instability of the vector field. Ignoring the spatial 
gradient term, the solution for the comoving electric field strength E evaluated on this 
background is given by 

where Eq is an integration constant, a is the scale factor and H is the Hubble expansion 
rate. When A is neg ative and -2XH 2 > M^, the energy density of the electric field 
can rapidly increase and eventually diverge, even when the expansion of the universe is 
accelerated. Our first goal is to take into account the back reaction of this growing vector 
field and examine the fate of the inflationary universe. 
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3 Dynamics of electric fields in axisymmetric Bianchi type I universes 

In this section and the next, we set M% = 1. 

3.1 Electric fields in an anisotropic universe 

In order to answer the question of back reaction, one needs to look at a fully non-linear 
system and solve both (2.5) and (2.7). The simplest generalisation of the FLRW universe 
that can accommodate a vector field is the axisymmetric Bianchi type I metric given by 

ds 2 = _ dt 2 + e 2a(t) [ e ~m) dx 2 + e mt) ( dy 2 + ^2) j ( g 

This metric is spatially flat. Note that the (mean) Hubble and shear expansion rates are 
given by 

H = a, a = $. 

We consider a homogeneous electric field along the x-direction, which in the gauge Aq = 
corresponds to the following coordinate basis components for the vector potential: 

A M = (0,A(t),0,0). 

The electric field strength seen by an observer moving with four- velocity = (1,0, 0, 0) is 
given by 

E(t) = -Ae- a+2 ?, 

where dots denote derivatives with respect to the comoving proper time t. We also include 
a perfect fluid with the equation of state p = (7 — l)p and constant 7. Hence, (2.5) yields 
the following: 



H 2 -a 2 



H + H' 



\(H + afE 2 + \E 2 + \p, 
b 3 



+- 



-2a 
XE 2 



l -E 2 - i(37 - 2)p -±(H + a)(H + 7a) E 2 



(3.2) 
(3.3) 



G 



A (H + a) 2 E 2 + ^E 2 - (7 - l)p + 4 (H + a) 



1 



1 



A 



a = -3Ha + -E z + - (H + a) (H + 7a) E 



XE 2 
6 



6 



A (H + a) 2 E 2 + ^-E l - (7 - 1) p + 4 (H + a) - 



1 



(3.4) 



We have already reshuffled the Einstein equations to put them into a convenient form. 
Eqs.(3.2) and (3.3) correspond to the Friedmann and Raychaudhuri equations, respectively. 
Eq.(2.7) can be written as 

2{H + a) 



E = — 



l + 2X(H + a) 2 
The fluid obeys the usual adiabatic decay law: 

p = — 37a/?. 



1 + 2X(H + a) z + 2X(H + a 



E. 



(3.5) 



(3.6) 
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Equations (3.2) - (3.6) form a closed set of non-linear ordinary differential equations. The 
Friedmann equation (3.2) measures the dynamical significance of each matter component. 

To gain an insight into the effect of the Horndeski's extra non-minimal coupling, let 
us assume that Horndeski's modification term is dominant, that is 



A (H + a) 2 E 2 



^E 2 ,p. 



The Friedmann equation (3.2) becomes 

H 2 -a 2 ~ X(H + <t) 2 E 2 , (3.7) 

which means that the universe must be strongly anisotropic when A < 0. Using the same 
approximation, (3.5) reduces to 

E 1 - 2Xa (H + a) 
E ~ A (H + a) ' 

Now, (3.3) becomes 

H + H 2 2a 2 , 

which is equivalent to the usual Einstein equation in an empty Bianchi I universe dominated 
by the shear because the back-reaction of the Horndeski term exactly cancels out at leading 
order. The same is true in the shear evolution equation, as (3.4) yields 

a ~ —3Ha, 

so that the universe isotropises in the same way as a (flat) universe containing only perfect 
fluids. We conclude that the Horndeski modification Ch is fairly innocuous despite the 
formidable appearance of its energy-momentum tensor. It should only be able to affect the 
evolution of the universe when its contribution is comparable to the conventional Maxwell 
term or the matter. 

3.1.1 Occurrence of a finite-time singularity 

The evolution equation for electric field (3.5) is rather similar to the linearized equation 
yielding the solution (2.8), although it is fully non-linear in the present setup. In fact, we 
can integrate (3.5) analytically and obtain 

This is essentially the same as the solution in FLRW background (2.8). Unless there is a 
mechanism within (3.2) - (3.6) that prevents (H + <r) 2 from reaching — 1/2A, there should 
be a range of initial conditions for which the system eventually hits a singularity. Since the 
system comes close to this singularity precisely when Horndeski's modified term becomes 
comparable to the Maxwell term 



A (H + a) 2 E 2 ~ E 



2 
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we may see unusual dynamical behaviours in this regime. 

Let us see what happens to the evolution of the spatial geometry when the system 
approaches the singularity. For this purpose, it is useful to write down the evolution 
equation for H + a, which we can cast into the following form; 



H + & 



3(F + ( r) 2 + (7-l)p 



+ 



1 + 2A (H + a) 



-E A 



(3.9) 



It is immediately clear that the right-hand side is negative definite when 7 > 1 and A < 
for initial conditions satisfying 



("+.) 2 >-i 



(3.10) 



Therefore, the singularity is inevitable regardless of the initial conditions for 7 > 1 as long as 
H + a > is sufficiently large initially. Although it does not apply to inflationary universes 
with < 7 <C 1, we already know H + a decreases monotonically while p dominates the 
evolution of the universe. Thus, if the condition (3.10) holds initially, E eventually grows 
and any matter domination, and hence inflation, comes to an end. Once the electric field 
starts to dominate the dynamics, the right-hand side of (3.9) is again negative definite and 
the singularity must be reached. While we are unable to eliminate the possibility that 
H + a turns to positive and the universe manages to avoid the singularity during a brief 
period of matter-electric equality, numerical calculations suggest otherwise (see figure 1). 




Figure 1. The approach towards the singularity for A = —2,7 = 0.1 and H + a > 0. The initial 
conditions arc p = 0.7, H = 0.4, a = 0.2. The universe is initially dominated by the matter with 
p/p = —0.9. There is no sign of avoiding the singularity located around t ~ 0.38. H becomes 
negative just before the singularity, which means the universe recollapses. 

Since the instability condition (3.10) roughly corresponds to the one for Horndeski's 
term to have a significant effect compared to the Maxwell term in the Lagrangian, it 
effectively rules out any sensible cosmological application of the theory with a negative A. 
We have not discussed the case of negative H + a since it describes either a collapsing 
universe or excessively anisotropic one. Close examination of (3.3) indicates recollapse or 
bounce right before hitting the singularity, both deriving from the violation of the weak 
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energy condition. These behaviours are also observed in our numerical solutions of the 
equations (see figure 2). 




Figure 2. The approach towards the singularity for A = —2, 7 = 0.1 and H + a < 0. The initial 
conditions are p = 0.7, H = —0.29, a = —0.2. The universe is initially contracting, but eventually 
bounces and rapidly expands before it reaches the singularity around t ~ 0.058. 



3.2 Expansion-normalised autonomous system 

In the previous subsection, we saw that the theory with a negative coupling, A, should 
lead to pathological dynamics when the Horndeski modification has an appreciable effect. 
On the other hand, when A is positive there is no danger of a singularity. The energy 
density is positive definite and we expect some viable cosmological dynamics. In order to 
carry out a more systematic investigation, it is always useful to rewrite the equations in 
terms of density parameters defined for each matter component, including the Horndeski 
contribution. It also enables us to apply the conventional methods of dynamical systems 
analysis. 

We introduce the following normalised variables: 

S = — , SIm = ^jj2> = 6A (H + a) 2 U M , ^ m = 3^2' 

where Qh correspond to A ^ respectively. The normalised Friedmann equation 

1 = s 2 + n H + n M + fim, (3Ti) 

will be used as the standard measure of the dynamical significance of each component. In 
particular, when > 0, all the parameters are bounded by 1 so that their values have 
a clear physical interpretation. Note that = pi/3H 2 for i = (M,H), where pu arid 
pH represent the minimal and non-minimal contribution to the vector field energy density, 
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respectively. From (3.3), we define the deceleration parameter as 



H 



+ 



2 (1 + S) (3fi A/ + fijy) 



3(7 - i)n m 
2(1 + e) z v 3ft M + n H 

(3 (3 + 5S) n w - (1 - E) «h) ■ 



(3.12) 



Following the standard method [78], we switch the time coordinate from t to a. Using 
(3.11) and (3.12), we derive the following evolution equations for the expansion-normalised 
variables: 



cZE 
da 
(Mm 
da 

dQ,H 
da 



d^l m 
da 



(q - 2) E - 2E 2 + n M ~ ^-^n m + 



2Q 



M 



q+l 



9. 



H 



1 + E 



M 



n 



H 



3(7 - 1R 

3Qm + &h J 3f2jv/ + VLh 



(1 



20 



ff 



9-2E + 



0// 6f?M — 0# 



n, 



1 + E 30 M + Oh 
3 7 -2 3( 7 -l)n H 



(1 + E) 



1 + S V 2 30 M + Off 

(2q -3 7 + 2) O m . 



(3.13) 
(3.14) 

(3.15) 
(3.16) 



These four equations are not independent since they are related by the first integral (3.11). 
3.2.1 Fixed points in the dynamical system 

We first classify the fixed points. Since the subsystem specified by 0# = is identical to 
the magnetic Bianchi type I discussed in [75, 76], we know there must be at least four fixed 
points of physical interest: 

Flat Friedmann universe : F 

(x,n M ,n H ,n m ) = (0,0,0,1). 



Electric Bianchi type I : E 

(E, Om, 0#, O m ) 



*r_i > ^(2- 7 ) (3 7 -4), 0,|(4 -7) 



The existence condition is 7 > 4/3. 
Kasner solutions : K± 

(E,0 M ,Otf,O m ) = (±1,0,0,0). 

In addition, there appears a fixed point describing a universe dominated by the Horn- 
deski energy density: 
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f2#-dominated universe : He 



{n,n M ,n H ,n m ) = (0,0,1,0). 

The deceleration parameter for this solution is q = 0, which is consistent with the 
analysis in section 3.1 where we showed that the back-reaction of the vector field 
exactly cancels out to leading order when the energy density is dominated by Horn- 
deski's modification term. 

There are subtleties regarding these fixed points. First of all, E and He do not 
represent physical space-times when A 7^ since they imply H = a = and (H + cr)" 1 = 
respectively. It does not mean these fixed points are irrelevant in the dynamics, however, 
since they may be reached asymptotically from finite H and a in the far past or future. 
We shall see an example of this in figure 7. Secondly, F and K± must be treated with care 
since some of the denominators appearing in the evolution equations (3.13)-(3.16) vanish on 
those fixed points. Again, it does not mean they are unphysical since they are well-behaved 
when appropriate limits are taken for the numerators. But the analysis requires evaluation 
of 0/0, which implies the stabilities may depend on the way the fixed point is approached. 
It is a consequence of the fact that H is not decoupled from the normalised variables, 
and implicitly appears in the definition of Qh- While a fixed point in the expansion- 
normalised equations usually represents a self-similar solution that is invariant under a 
scale transformation, the dynamical effect of the Horndeski modification depends on the 
scale of the curvature, or the size of the universe. Therefore, it is not surprising to see the 
stability change depending on each orbit with its specific value of H. We will find this is 
indeed the case. 

3.3 Dynamics around the matter-dominated solution 

From physical point of view, by far the most interesting fixed point is F because it can be 
regarded as a model of the late-time evolution for the universe when 7=1 (dust) or 7 = 4/3 
(radiation), and also a model of inflation when 7 < 2/3. We have already mentioned the 
instability against perturbations of the electric field for A < 0. The condition for the 
occurrence of a singularity (3.10) translates into 

3n M + ^H<0 (3.17) 

in the new variables. Here, we shall see that this condition coincides with the instability 
condition for F and otherwise the dynamics is trivial. We also study the stability for 
positive A and show that it depends on the value of H. 
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3.3.1 Linearisation 

For the purpose of linearisation around F, it turns out to be convenient to eliminate Q m 
using the Friedmann equation (3.11) and rewrite the equations as 

— = Ms,n M ,n H ), (3.18) 

da 

f M (x,n M ,n H ), (3.i9) 



da 

dn H 

da 



f H (x,n M ,n H ), (3.20) 



whose right-hand sides we avoid writing down explicitly as they are lengthy. While there 
are apparent 0/0s in those equations when they are evaluated on F, they should be all well 
defined if an appropriate limit is taken along an arbitrary reference orbit. To proceed, we 
evaluate the functions fi,i = S, M, H for S = and then take the limit (fi^, Q#) — > (0, 0): 

lim s h(0,n M ,n H ) = -2~fL 2 , 

lim / m (0,1]m,^) = -27L2, 

lim f H (0,ttM,tt H ) = 67L2, 
(fi M ,n H )->(o,o) 

where we have introduced a notation 



Ln = lim 



(n M ,n H )-Ko,o) 3Q M + &h ' 

which will also be used later. When Q,h > 0, or equivalently A > 0, we have L2 = 
and the fixed point F is always well defined as it should be. The existence of the limit is 
inconclusive when < (A < 0) and the orbit satisfies 

lim = -3. (3.21) 

(n Ml n H )-Ko,o) fijw 

However, such an orbit merely represents one that ends up in the singularity 2\(H + cr) 2 = 
— 1. Since we have already discussed this case in detail, we exclude those orbits from our 
consideration here. As long as an orbit does not hit the singularity when approaching F, 
L2 should exist and be equal to zero. 

Since the right-hand sides can be evaluated only as a limit associated with each ref- 
erence orbit, the linearisation takes an extra step. We first expand the equations around 
an arbitrary point (So, &mo, Qho) and then take the limit (So, &mo, &ho) (0, 0, 0). The 
resultant linear equations are given as follows: 

— = -(7 - 2)5S + (2 + 67L 2 ) 6U M + - [7 + 2 + 4 7 Lx (L x - 2)] 8n H , (3.22) 

(37 - 4 + 6 7 L?) sn M + 27 (Li - 1) 2 sn H , (3.23) 



da 

-18jL'(Sn M - 2 [2 + 37L1 (Li - 2)] 5n H , (3.24) 



da 
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where the 6s preceding the variables denote their small perturbation. Unless the orbit is 
the singular one specified by (3.21), L\ is finite and therefore these linearised equations are 
well defined. However, the asymptotic value of L\ does depend on each orbit. Going back 
to its definition, one notices that 

u n 



1Z + 3 
where 

I'M 

We already know its behaviour near F since we have 

n ~ 6XH 2 ~ 2\p. 

Solving (3.6), we obtain 

U = K e- 3 ~ (a , (3.25) 
where IZq is an orbit-specific constant. We also note that 

6n H = K6n M (3.26) 

along each orbit, so that we have an additional linear constraint. Now, the linearized 
equations (3. 22)- (3. 24) can be written 

l( 1 -2)6E + j^\6 + ^n[l0 + 3 1 -n(3j-2)]\6n M , (3.27) 



da 2 

6-fU 



d6VL M 

— = dllM 

da 

d6Q,H 



da 



6n 



H 



3 7 -4 + 
-4 + 



3 + 11 
6jK 



3 + 11 



(3.28) 
(3.29) 



and we can easily read off the eigenvalues of the linearization matrix 



^( 7 -2),3 7 -4 + J^-,-4 + J^-V (3.30) 
2 W ' 1Z + 3' 11 + 3) v 1 

Notice that the eigenvalues are orbit and time dependent through 1Z(a). As it was necessary 
to take a non-standard approach to obtain the eigenvalues, we will later confirm the validity 
of the result by performing numerical calculations. 

3.3.2 Stability of the matter-dominated solution 

The first of the eigenvalues (3.30) is negative and represents the stability of F against per- 
turbation of S. Since the third eigenvalue is always smaller than the second, the condition 
for the stability is 

Let us first consider A > 0, which corresponds to 1Z > 0. In this case, we have 

1Z 

0< nTs <1 ' 
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and consequently F is definitely stable for 7 < 4/9, which means positive A cannot be 
relevant in the context of inflation, and F is unstable for 7 > 4/3. For 4/9 < 7 < 4/3, the 
stability is orbit-dependent. When an orbit satisfies 

it runs away from F. When TZ is smaller than this threshold value, the orbit is attracted 
towards F. Note that the value of TZ is time-dependent so that the stability can change 
over the course of the evolution. In particular, from (3.25), TZ is monotonically decreasing 
as long as the orbit stays close to the matter-dominated solution. We can immediately 
conclude that for orbits with TZ < (12 — 97) /(97 — 4), the stability does not change as the 
universe expands. For those satisfying (3.31), they typically become stable asymptotically 
in the future since TZ can only increase when the universe is dominated by both Maxwell's 
and Horndeski's terms. For a physically interesting range of initial conditions, we shall 
later confirm that the instability of the electric field saturates before the orbit goes too far 
from F and eventually comes back to it. 

For A < 0, the dynamics is very different - depending on 1Z —3. As was already 
mentioned, the critical value 1Z = —3 corresponds to the singularity 2\{H + a) 2 = — 1 
discussed in section 3.1.1. Firstly, 1Z G (—3,0) implies 72-/(72. + 3) £ (— 00, 0) and therefore 
the orbits in this range are stable as long as 7 < 4/3. For 1Z < —3, we have TZ/ (TZ + 3) G 
(l,oo) and TZ monotonically increases in the vicinity of F. It approaches TZ = —3 from 
below and therefore any orbit eventually becomes unstable. Since we already know TZ = —3 
is the singularity, we conclude that the orbits with this range of initial TZ can never settle 
down at F, regardless of the equation of state parameter 7; see figure 3 for simulation 
of an inflationary universe with initial conditions satisfying TZ = —100. The solution 
approaches F until TZ is close to the critical value in which case the universe moves away 
from F. After leaving F, we expect the system enters the regime where f2 m is dynamically 
negligible. Then the analysis in the section 3.1.1 applies and the singularity is inevitable. 

3.3.3 Dynamics with A > in a radiation- or dust-dominated universe 

We have found that the theory is quite innocuous for positive non-minimal coupling con- 
stant, A, in which case there is no singularity and the Friedmann solution F is stable at 
late times. However, as shown above, there is a transient period where F is unstable for 
fluids satisfying 7 > 4/9. In this period, Qm and £Ih grow and, if the instability does 
not saturate before they become too large, the universe will eventually become strongly 
anisotropic. Since this introduces potential problems in the radiation or dust-dominated 
epoch, it is of interest to specify the range of initial conditions such that the universe is 
close to F at all the subsequent times, i.e., \£l m — 1| <C 1. In this subsection, therefore, we 
investigate the dynamics close to F in more detail for dust (7 = 1) and radiation (7 = 4/3) 
with a positive non-minimal coupling constant (A > 44> VLh > 0). 
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Figure 3. Simulation of (3.13)-(3.15) for an inflationary universe with 7 = O.f and initial conditions 
close to F (S = -O.Of , tt M = 0.0001, O ff = -O.Of). The singularity at 3Q M + = occurs at 
time a — 1 1 .65. 

We can easily solve the linearized equations (3.28) and (3.29) exactly: 

sn M (a) = sn M0 ^|±^V e -(*-*y)«, (3.32) 

6Q H (a) = Sn H0 (|^) 2 e" 4a , (3.33) 

where 5Qmo and 5Qho are integration constants. At late times, when 1Z oc e _37 ° — > 0, 
the Maxwell and Horndeski densities decay as SQm oc e~^~ 3 ^ a and 5£1h oc e~ Aa . This 
is consistent with F being an attractor at late times as shown above. When 1Z S> 3, the 
Maxwell and Horndeski densities grow as 5Qm oc e^ 7-4 ^" and 5Qh oc e 2 ( 37 ~ 2 ) Q . We note 
that the maximum values of 5VLm and 5£Ih occur when 1Z = (12 — 9j)/{— 4 + 97) and 
1Z = 6/(37 ~~ 2), respectively. Taking into account the linear constraint (3.26), it follows 
that the maximum values of both Qm and £Ih are roughly equal. Now we find that the 
initial conditions must satisfy 

Sn m ^ (SnMof 2 O X^iEoHo)- 1 (3.34) 

in a radiation-dominated universe and 

<K2jto < (anjifo)* & X <. (E )-HH )-I (3.35) 

in a dust dominated universe to ensure that |fi m — 1| <C 1 at all the subsequent times. In 
figures 4 and 5, we show a numerical integration of the full non-linear equations (3.13)- 
(3.15) for radiation and dust, respectively. Since the initial conditions just barely satisfy 
the conditions (3.34)-(3.35), the peak values of Qh and Qm are at the one-percent level 



- 14 - 



of the total energy budget. Note that, to good accuracy, we have £ when 5s 
fin- This introduces the possibility of cancelling the effects of spatial anisotropy on the 
Cosmic Microwave Background (CMB), which will be discussed in section 5. Under the 
assumption that the theory describes a generalised electrodynamics, we show in section 5 
that the amplification of the electric field must come to an end before the start of big bang 
nucleosynthesis, when the temperature is T ~ 1 MeV. This still leaves the possibility open 
for a huge amplification of electric fields in the period between inflationary reheating and 
the nucleosynthesis. In the period of amplification, the electric field grows very quickly, 
Q.M oc e 8a . After the peak value is reached, Qh rapidly decays and soon becomes negligible. 
At that stage, the dynamics becomes similar to the conventional electrodynamics; Qm 
decays logarithmically (constant at the linear level) until the dust-dominated epoch when 
it decays as Qm oc e~ a . 





Figure 4. Simulation of (3. 13)- (3. 15) for radiation (7 = 4/3) and a positive non-minimal coupling 
coupling constant (A > 0). Initial conditions (£ = —3* 10 -9 , Qjvf = 10~ 14 , fin = 2* 10 -8 ) are such 
that the orbit is always close to the Friedmann solution F. From the logarithmic plot to the right 
it is clear that the ratio Qh/^m is monotonically decaying in agreement with equation (3.25). 



3.4 Stability of the other fixed points 

Given the complexity of the dynamics, it is also helpful to analyse the stability of the other 
fixed points. Here we present the eigenvalues for K + , M and He- The linear stability 
analysis for K- is inconclusive since we have H + a = 0, E = there, and so the Horndeski 
energy density 3X(H + a) 2 E 2 is generically second order in perturbation. Our numerical 
simulations in the next subsection indicate that K— is a past attractor for A > 0. 

When either Qm or Qjj is nonzero, the stability analysis is straightforward since there 
is no orbit-dependence. We obtain the following eigenvalues: 



Fixed point M 



2 - 7 ± y/(2 - 7 )(3 7 2 - 177 + 18) 



-3 7 
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Figure 5. Simulation of (3.13)-(3.15) for dust (7 = 1) and a positive non-minimal coupling coupling 
constant (A > 0). Initial conditions (S = —7.5 * 10 -8 , Qm = 10~ 14 , = 5 * 10~ 7 ) arc such that 
the orbit is always close to the Friedmann solution F. From the logarithmic plot to the right it is 
clear that the ratio VLh/^m is monotonically decaying in agreement with equation (3.25). 



Fixed point He 

(±2,2-37). 

Whenever M exists (7 > 4/3), it is a future attractor. The new fixed point He is a 
saddle regardless of 7, and the orbit is temporarily attracted towards it for initial conditions 
close to F that fail to satisfy the conditions (3.34)-(3.35). We did not find the dynamics 
around this solution to be of any phenomenological interest. 

The eigenvalues of K + is dependent on each reference orbit and the linearisation can 
be carried out in a similar way as for F in section 3.1.1. In the end, we derive the following 
equations: 

^ _ 3(2 - 7M £ + 24-^ + 7^ + 38-3^ + 4)] 
da y U 2(ft + 3) 

dsn M 10ft - 6 



da ft + 3 

d5n H 4K-24 



da TZ + 3 

and read off the eigenvalues: 

10ft -6 4ft -24 



3(2-7), 



ft+3 ft+3 



The first eigenvalue is associated with the perturbation of VL m and positive so that K + 
cannot be a future attractor. The time-dependence of ft = 6X(H + a) 2 around K + is given 
by evaluating (3.9) with H ~ a as 



ft ~ ftoe 



-6a 
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For 1Z < 0, the situation is analogous to F except that K + is irrelevant to future asymptotic 
behaviour. An orbit with 1Z < —3 initially reaches the singularity while K + is a saddle 
point for the others. For 7Z > 0, the stability changes over the course of the evolution. 

3.5 Numerical Analysis 

We conclude this section by showing two visualizations that confirm the analysis of this 
section, and suggest that K- is a past attractor (which is the case in the minimally coupled 
theory) for A > 0. Figure 6 shows a phase portrait of the invariant subset Q, m = (a 
universe without the fluid). This subset is effectively two dimensional by the Friedmann 
equation (3.11) and the phase space is completely characterised by the variables £ and £Im- 
The green shaded region corresponds to A > 0. It clearly shows that K_ is a past attractor 
for A > on this subspace. Note the existence of stream lines directed towards in 
the A < region indicating orbit dependence of the stability for a negative non-minimal 
coupling constant. In the region A < 0, there are two distinct flows inside and outside 
the bold red lines that denote the location of the singularity. The outer region is marked 
with a red mesh and corresponds to the pathological region (3.17) where the singularity is 
inevitable. Note that (3.17) can be written 1 + 2Qm < £ 2 + which explains why the 
entire region £ 2 < 1 is non-pathological on the subspace Q, m = 0. When a fluid is included, 
however, the singularity can be reached even from a small initial shear (for example figure 
3). Figure 7 shows the full three-dimensional phase flow for a few orbits with A > 0, which 
demonstrates the past-stability of K_ and confirms the time-dependent stability of F. 

4 Dynamics of a magnetic field in axisymmetric Bianchi type I 

In the usual electromagnetism, Maxwell's equations treat electric and magnetic fields in 
a symmetric manner. In the context of Bianchi cosmologies, source-free pure electric and 
magnetic fields are mathematically indistinguishable. The modification Ch breaks this 
duality. Let us consider a homogeneous magnetic field along x-axis in the spacetime (3.1). 
In parallel to the electric field of the previous section, we define B as the magnetic field 
seen by a comoving observer: 



Contrary to the electric case, equation (2.7) is trivially satisfied for this Faraday tensor. 
Instead the evolution of the comoving magnetic field is given by the Bianchi identity dF = 
which leads to 



In contrast to the peculiar dynamics of the electric field, there is no modification to the 
evolution equation for the magnetic field since it comes from the closedness of the field- 
strength 2-form. Then, we can solve it easily to obtain 




(4.1) 



B = —2 {H + a) B. 



(4.2) 



B = 



Bo 



(4.3) 



e 2(<*+/3) • 
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Figure 6. Phase flow for the subsystem il m = in the electric case. The green shaded region 
corresponds to A > 0. The red bold line is the position of the singularity, while the red mesh is the 
region where the singularity is inevitable. 



This simply means the magnetic field is adiabatically decaying due to the expansion of the 
universe. 

The Einstein equations can be written in the following convenient from: 



H 2 



H + H z 



1 



1 



a ~\~ ~p ~\~ -B - -X (H + a) (H - 2a) B\ 



+ 



1 + XB 2 
1 



H 2 + Ha + 2a 2 + - (7 - 1) p - 2XH (H + 2a) B 2 
3 

7-1 



H + a 



2(1 + \B 2 ) 2 
1 



(H + af + l-^ P - l -B 2 -\x(H + afB 2 

6 6 



1 + XB 2 



3a (H + a)- 3XH (H + a) B z + 



(4.4) 
(4.5) 



(4.6) 



Again, there appears to be a problem for negative A when XB 2 ~ — 1. This time, the 
singularity stems from the Einstein equations instead of the evolution equation for F a b as 
in the electric case. From (4.6), it is clear that H + a is positive definite in an expanding 
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Figure 7. Simulations of the non- linear equations (3. 13)- (3. 15) for a radiation fluid (7 = 4/3). 
All the orbits correspond to a positive non-minimal coupling, A > 0. The thick red line starts from 
K + , gets repelled from F and eventually comes back to it, which confirms the time-dependence of 
the stability for F. All the other orbits traverse from K- to F indicating K- is a past attractor. 

universe H > for initial conditions satisfying 

H + a>0, -XB 2 > 1. (4.7) 

According to (4.3), B 2 decays in this regime and consequently the system has to reach 
the singularity XB 2 = — 1. As one can see in figure 8, the shear becomes negative before 
reaching the singularity. The behaviour is insensitive to the initial conditions or value of 
A whenever —XB 2 > 1 initially. This condition coincides with the one for the Horndeski 
modified term to have a significant contribution to the dynamics. We conclude that A < 
is pathological for the magnetic case too. 

While there is no divergent behaviour for positive A, the energy density of the magnetic 
field is not necessarily positive, in contrast to the electric field case. This will not cause 
any problem for ordinary expanding universes with H + a > since the solution (4.3) 
ensures monotonic decrease of B. However, it may result in unusual behaviours when we 
go backwards in time. We repeat the dynamical system analysis of the previous section 
and show that A > is cosmologically viable. 

4.1 Expansion- normalised autonomous system 

We introduce the following normalised variables: 

S = -|, n M = -^, n H = -^X(l + E)(l-2E)B 2 , 
The normalised Friedmann equation takes the canonical form 

i = z 2 + n H + n M + n m . (4.8) 



n, 



m 2 ' 
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Figure 8. Occurrence of the singularity for A = —2,7 = 1. Initial conditions arc H = 5.1, a = 
2.9, B = 1. In fact, the dynamics is more or less the same for any values of parameters or initial 
conditions as long as —\B 2 > 1. In this case, while the magnetic field is subdominant in the 
Friedmann equation (4.4), it still affects the evolution of H and a and causes their divergences. 
Note that a/H — > — 1 as the singularity is approached. This implies the approach towards the 
singularity appears as a flow into in the expansion-normalised variables. 



In contrast to the electric field, is not positive definite regardless of the sign of A. 
Following the steps of the previous section, we rewrite (4.5) as a defining equation for the 
deceleration parameter: 



2 (1 + S) (1 - 2E) 
'3^-2(l + E)(l-2E) 

2(1 + E) 2 (1-2E) 2 
[Sn H - 2 (1 + S) (1 - 2E)] 5 



l + S + 2S 2 + 2( 7 -l)O m + 



3(1 + 2E) 



(l + E) 2 + ( 7 -lR 



(1 + E)(1-2E 
2(1 + E) 



-n 



M 



1 - 2E 







H 



We derive the following evolution equations for the normalised variables: 



dT, 
da 



(1 + <7)(1 + E) 

6(1 + E)(1-2E) 



+ 



M 



da 
diln 
da 

dQrn 

da 



3n H — 2 (1 + E) (1 — 2E) 
2(<7-l-2E)fi M , 

1 + 4E 



S(1 + S) + I^ + 27T 1 2E) 



n 



-4(i + z)n H - 

(2q -37 + 2) n r 



(1 + E)(1-2E) H da' 



4.2 Fixed points and their stabilities 



(4.9) 

(4.10) 
(4.11) 

(4.12) 



We repeat the standard stability analysis. The structure of the state space is analogous to 
the electric case. There are four fixed points residing in the conventional magnetic Bianchi 
type I, and another with non- vanishing Qh- 
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Flat Friedmann universe : F 



(E,n M ,n H ,n m ) = (0,0,0,1) 

with eigenvalues 

(_4,-|(2- 7 ),37-4). (4.13) 

Since the evolution equation of the magnetic field is well-behaved, the linearisation 
can be carried out as usual. The stability does not change from the usual magnetic 
cosmologies. F is a future attractor for 7 < 4/3. Note the zero eigenvalue in the 
radiation case (7 = 4/3). 

Magnetic Bianchi type I : B 

(E, M , Sl H , n m ) = (j 1 ^, ^ (2 - 7) (37 - 4) , 0, 1 (4 - 7) 
with eigenvalues 

(~ (2 - 7 ± V(2 - 7)(37 2 - 177 + 18)) , -3 7 ) . (4.14) 

The existence condition is 7 > 4/3 and it is always stable. 
Kasner solutions : K± 

(z,n M ,n H ,n m ) = (±1,0,0,0). 

The eigenvalues for K + are easily computed as 

(3(2 -7), -2, -8), (4.15) 

indicating it is a saddle while being a future attractor for the subsystem f2 m = 0. 
For K_ , we have the same indeterminacy of the linear stability as was encountered 
in the electric case. It will be examined by numerical analysis later. 

fijf-dominated universe : Hb 

(£, Qm, ^m) = ^' 4' ^ 

This fixed point does not itself represent a physical solution since it corresponds to 
the limit —X(H + cr)(H — 2a) — > 00. As in the electric case, this does not mean 
that it is dynamically irrelevant since it can be reached asymptotically from finite 
H and a (see figure 11). Note that q = for this fixed point which means that the 
back-reaction of the vector field exactly cancels out like for the Horndeski dominated 
solution in the electric case. Its eigenvalues are given by 

(-12, -3(7 + 2), 6). (4.16) 

Therefore it is always a saddle point. 
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4.3 Numerical analysis 

Since the uniform magnetic field decays adiabatically, the evolution close to the Friedmann 
solution F is trivial compared to the electric case. Near F, £Im evolves as in the con- 
ventional electromagnetism; logarithmically decaying (constant at the linear level) in the 
radiation era (7 = 4/3) because of the zero eigenvalue, and in proportion to e~ a in the dust 
era (7 = 1), see [79-81]. It follows immediately from its definition that oc B 2 ex e -4 ", 
i.e., the evolution is independent of the equation of state of the perfect fluid. See figure 9 
for a simulation of the full non-linear equations (4.9)-(4.11) close to F. 

It is again helpful to visualise the global structure of the phase space. Figure 10 shows 
a phase portrait of the invariant subset Q OT = 0. The green shaded region corresponds to 
A > 0. It is clearly seen that K_ is a past attractor for A > on this subspace. Note 
the existence of stream lines directed towards K_ in the region A < 0. They indicate orbit 
dependence of the stability for a negative non-minimal coupling constant value. The red 
mesh corresponds to the pathological region (4.7) where a singularity, marked with the red 
bold line, is inevitable. Note, from its definition, that Qh = when E = 1/2. Therefore, 
although it has no physical significance, the point (fijvr,E) = (3/4,1/2) appears as if it 
were a fixed point. Figure 11 shows the full three-dimensional phase flow for a few orbits 
with A > 0. Overall, our simulations indicate that K_ is a past attractor for A > 0. 



l r 1 r 




468 10 2468 10 



Figure 9. Simulation of (4.9)-(4.11) close to the Friedmann solution F for dust (7 = 1) to the left 
and radiation (7 = 4/3) to the right. Initial conditions: S = 10~ 8 , VLm = 10~ 8 , = 10~ 2 . 

5 Observational constraints 

In this section, we shall discuss the possibility of detecting any signature of the theory 
with a positive coupling constant A in the standard cosmology. Since the correction term 
comes with a factor of the space-time curvature which is tiny in the units of Planck mass, 
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Figure 10. Phase flow for the subsystem f2 TO = in the magnetic case. The green shaded region 
corresponds to A > 0. The red bold line is the position of the singularity, while the red mesh is the 
region where the singularity is inevitable. In the meshed region, the flow goes towards K-. For 
A > 0, K_ appears to be a past attractor. 

one expects that a large value of the parameter A is required to see any appreciable effect. 
Combining the constraint on spatial anisotropy coming from big bang nucleosynthesis [82] 
and the validity of the Maxwellian electromagnetism in terrestrial environments, we will 
show that the growth of the electric effect discussed in section 3.3.3 must have come to an 
end before the CMB temperature dropped below T ~ 0.1 MeV. After that time, VLh decays 
much faster than Q^Mi an d so the effect of the Horndeski Lagrangian term is completely 
negligible at late times. Since we restrict ourselves to the dynamics around FLRW in this 
context, the magnetic case will not be discussed in detail due to its trivial evolution seen 
in the previous section. 

A non-minimal coupling of electrodynamics to gravity is known to cause a number of 
observational effects, such as frequency-dependent bending of light in gravitational fields 
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Figure 11. Simulations of the non-linear equations (4.9)- (4. 11) for a radiation fluid (7 = 4/3). 
All orbits correspond to a positive non-minimal coupling, A > 0. For each orbit, phase flow is in 
the direction towards the Friedmann solution F, and initial conditions are close to K_ with non- 
vanishing perfect fluid and the vector field, which underwrites the past stability of K- . The thick 
blue orbit temporarily leaves the boundary of the box. After it returns, it has transient periods 
close to Hb and K + before settling down at F. 



. Since those effects must be suppressed in terrestrial environments, we require 



\C H \ « 



(5.1) 



evaluated around the Earth. In an orthonormal basis, the components of the Riemann ten- 
sor for the Schwarzschild solution are of order R a b c d ~ Rs/r 3 \ where R s is the Schwarzschild 
radius and r is the radial coordinate in Schwarzschild coordinates. Inserting the radius and 
mass of the Earth, we find that the constraint (5.1) is equivalent to 



|A|« 



R. 



■Ml ~ 10 90 , 



(5.2) 



in which the upper limit appears enormous since we have used a large mass scale (M p i) 
to normalise the non-minimal coupling constant (which is a reasonable convention in the 
context of modified gravity). From the linear analysis in section 3.3.3, we know that both 
Qm and Qh start to decay when 



A 



H 2 



(5.3) 



is close to unity. At the time 1Z ~ 1, the expansion rate must satisfy 



H > 10" 45 M 



pi 



(5.4) 
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to be consistent with the constraint (5.2). Using the approximation T ~ (HM p i)z, where 
T is the CMB temperature, this imposes the following condition on the temperature at 
K~ 1: 

T n ^i > 10~ 23 M pl ~ 0.1 MeV. (5.5) 

Thus, to be consistent with the local constraint on terrestrial electromagnetism, the growth 
of the normalised vector field energy is required to end at latest around the time of big bang 
nucleosynthesis when T^bn = 0.1-1 MeV. Since 1Z is monotonically decreasing in the perfect 
fluid dominated universe (TZ oc e _4a in the radiation-dominated epoch and TZ oc e~ 3a in the 
dust-dominated epoch), we set 1l< 1 for T <C 1 MeV. It follows that the energy density 
of the vector, both in the electric and the magnetic case, will evolve as in the conventional 
electrodynamics, i.e. E 2 ,B 2 oc e~ 4a and therefore Qh oc e~ 4a . Using that T oc e~ a , we 
obtain the relation 

^ bn -^(^f) 4 , (5-6) 

where and O^ 11 are the values of Qh today and at big bang nucleosynthesis, respec- 
tively. The CMB temperature today is T ~ 10~ 4 eV. We therefore have tt^ 11 ~ lO 40 ^. 
The theory of big bang nucleosynthesis requires that the universe was dominated by the 
CMB and three effectively massless neutrino species at Tbb n; i- e - that Q}^ n must be smaller 
than unity. This corresponds to the following constraint at the present time: 

\n° H \ < 10~ 40 . (5.7) 

We conclude that the effect of Qh on the evolution of the universe becomes increasingly less 
significant after the big bang nucleosynthesis. It is interesting to compare this constraint 
to that on a uniform magnetic field in the conventional electrodynamics. In that case the 
strongest limit on a uniform magnetic field today comes from the CMB anisotropy and is 
[84] 

B < 10" 9 Gauss (5.8) 

or equivalently 

n o = go < io- 12 , (5.9) 

which is much weaker than the upper bound on the Horndeski part of the Lagrangian (5.7). 
Evaluating these constraints at decoupling time (z ~ 1100) 

\n<%\ ~ n^noo 4 < io- 28 , 0$ ~ n° M noo < io~ 9 , (5.10) 

we note that <C flj^- Consequently the additional anisotropic stress induced by 

the non-minimal coupling is completely negligible in the entire period between decoupling 
time and today. Now, since the upper bound on a uniform magnetic field comes from the 
modified temperature pattern in the CMB created by the shear (which is sourced by the 
anisotropic stress) , it is clear that the constraint so far obtained is strong enough to exclude 
any signature of £Ih in the CMB, and (5.9) holds also in Horndeski's generalised theory. 
This result reflects the fact that decays much faster than Qm after reaching 1Z ~ 1. 
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In a similar manner, even a direct detection of homogenous cosmological magnetic fields 
would not improve the constraint on A either since the regime £Ih > must come prior 
to big bang nucleosynthesis. 

So far we have assumed that the vector field is identified with photons, i.e. that the 
model is an extension of the ordinary electrodynamics. If this is not the case, and the 
model is merely taken to be a hypothetical vector-tensor theory, the bound on the non- 
minimal coupling constant (5.2) is not valid and the normalised vector field energy can 
grow also at late times if A is enormous. To good approximation, we have £ ^ when 
Qm Ujj, as noted in section 3.3.3, and this could have some interesting phenomenological 
consequences. Since the AT/T created in the CMB by the shear is an integrated effect, 
which is proportional to [85] 



a cancellation effect will occur if the peak value of the normalised vector field energy 
(around 1Z ~ 1) occurs between decoupling time (adc) and today (ao). First, one will have 
a period with £Ih > where both Q,m and VLh grow and £ < 0, followed by a period 
with £Ih <^m where both Qm and Qh decay and £ > 0; see figure 5. By fine-tuning the 
non-minimal coupling constant, the cancellation can be made exact leaving no effect on the 
CMB despite having a relatively large amplitude of shear. Numerical simulations show that 
this cancellation requires that Q,h peaks at a low redshift; as an example, Qh must peak at 
redshift z ~ 0.6 for the initial conditions in figure 5. The shear still needs to be consistent 
with the supernovae la data which is isotropic to the 1% level, but this is significantly 
relaxed compared to the CMB bounds in models where £ is positive or negative definite 
where |£| is never greater than ~ 10 -5 . The possibility of such a cancellation effect was 
mentioned in [86, 87] and to the best of our knowledge, Horndeski's theory is the first 
instance of a Lagrangian for which such a scenario may be dynamically realised. 

6 Conclusion and outlook 

In this work, we have explored the cosmological consequences of a most general vector- 
tensor theory that gives second-order field equations equivalent to Maxwell's equations 
when evaluated in flat space-time. We have focused on a simple dynamical setup which is 
well understood in the conventional electromagnetism, namely the case of a uniform electric 
or magnetic field in an axisymmetric Bianchi type I universe. We have investigated the non- 
linear evolution equations by using several different approaches, including the conventional 
dynamical systems analysis in terms of Hubble normalised variables such as £Im and £Ih- 
In contrast to the conventional electromagnetism, the theory treats electric and magnetic 
fields in different ways and we have discussed each case separately. 

It has been found that the cases of a positive and negative non-minimal coupling pa- 
rameter A are drastically different. For A < 0, we identified physical finite-time singularities 
where the deceleration parameter diverges. In the electric case, the singularity comes from 
the modified evolution equation for the electric field, while in the magnetic case, it stems 
from the Einstein equations. In both cases, we identified a range of initial conditions (which 
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depend on A) for which the singularities are inevitable. Although we have been unable to 
show the singularity is unavoidable during inflation, its mere existence in state space is a 
problem in itself and suggests that the theory is pathological for A < 0, regardless of the 
other matter fields in the universe. Our result effectively rules out the possibility of gener- 
ating anisotropy or magnetic field during inflation by employing this type of non-minimal 
coupling. 

For A > 0, it appears that the modification is relatively harmless as there is no singular 
behaviour and the Friedmann solution (F) is the future attractor for any perfect fluid 
satisfying 7 < 4/3. For inflationary universes with 7 <C 1, we have found F to be stable, and 
consequently all anisotropics are washed out. In radiation and dust-dominated universes, 
however, we find a rather interesting phenomenology. Namely, in the electric case, the 
stability of F is time dependent and there is a transient period where both VLm and Qh 
grow. Both Om and reach their maximum values around the time Q,m ~ and then 
Q.M decays adiabatically but the decrease of is much faster. Demanding the validity of 
the Maxwellian electromagnetism in terrestrial environments, this peak value has to occur 
at the latest by the time of big bang nucleosynthesis, which means the cosmological effect 
of VLh afterwards is always negligible. Therefore, all the cosmological constraints after 
nucleosynthesis are based on the dynamics of £1 m and hence cannot be used to improve the 
upper bound on the non-minimal coupling constant A. Nevertheless, a huge amplification 
of electric field energy may have taken place in the early universe (&m °c c 8q when F is 
unstable). Even if electric fields are washed out during inflation (Qm °c e~ 4a if 7 ~ 0), 
there could be enough time between reheating and big bang nucleosynthesis for jy/ to grow 
significantly if inflation lasts not much more than 60 e-folds and the reheating occurs at a 
relatively high energy scale. After the amplification saturates, Qjj decays quickly and soon 
becomes negligible so that it does not cause any significant problems at later times. We 
note that this amplification is characteristic of uniform electric fields; for uniform magnetic 
fields, the dynamics near F are practically identical to the conventional electromagnetic 
case. 

Contrary to the minimally coupled Maxwell case, we found that the generalised energy- 
momentum tensor may source both a positive and negative shear. If the theory is merely 
taken to be some hypothetical vector-tensor theory (not electromagnetism), so that the 
Hubble normalised field energy may grow at late times, this opens up the possibility for 
cancellation effects in the CMB. In particular we saw that if f2# peaks around redshift 
z ~ 0.6, all CMB effects of the anisotropic expansion cancel exactly. The possibility of 
such a cancellation effect was mentioned in [86, 87] and to the best of our knowledge, 
Horndeski's theory is the first instance of a Lagrangian for which such a scenario may be 
dynamically realised. 

In this work we have neglected all spatial derivatives, and focused on understanding 
the full non-linear nature of the system. Physically, uniform electric and magnetic fields 
correspond to a long wavelength limit of cosmic electric or magnetic fields. This approach 
has been fruitful as it enabled us to identify singularities at the background level and 
therefore rule out any interesting cosmological application of this theory for a negative 
non-minimal coupling. Although this makes the theory irrelevant in the context of in- 
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flation, we have demonstrated that a non-trivial phenomenology is possible in the early 
radiation-dominated universe. We point out that the short-wavelength limit appears to be 
an interesting direction of further investigation. Especially, when spatial derivatives are 
taken into account, it would be interesting to see if an amplification of cosmic magnetic 
fields could be possible before big bang nucleosynthesis. It is also necessary to check if the 
model avoids instabilities from inhomogeneous perturbations. 
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